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Abstract. From medium-resolution radio images, DR 16 was suspected to be a large cometary nebula. To test 
this suggestion we obtained a higher resolution (15") VLA continuum map. We also analyzed data from the 
Canadian Galactic Plane Survey in continuum, Hi line, and IR. These data were supplemented by published 
near-infrared (J, H, K) stellar photometric results and MSX 8.28 )im data. We suggest that DR 16 is the diffuse 
Hn region of an ongoing star formation site at a distance of about 3 kpc. The complicated radio picture arises 
from the superposition of diffuse H II with the remains of a giant bipolar outflow. The outflow was generated by a 
probable Herbig AeBe star, and the lobes are the remnants of its working surfaces. Additional ring-like features 
are discussed. DR 16 is part of a larger volume of space in the local spiral arm where star formation is an ongoing 
process. 

Key words. Stars: early-type - stars: emission-line, Be - ISM: general - (ISM:) Hn regions - ISM: individual 
objects: DR 16 - Radio continuum: ISM 



1. Introduction 

In the direction of Cygnus X our line of sight is tangen- 
tial to the local spiral arm, and the superposition of many 
objects makes the region one of the most complex in the 
sky. Although the radio emission in Cyg X is mainly ther- 
mal in origin, individual sources reveal a surprisingly rich 
spectrum of properties, with objects in different stages of 
evolution occupying volumes of space with very different 
star formation histories. Uncovering the nature of individ- 
ual objects and establishing their distances demands ex- 
tensive and thorough observations. The most recent radio 
continuum surveys (Wendker et al. 1991 Paper XVIII) 
showed that a resolution of about 1' begins to resolve 
many sources, so that observations of higher resolution 
will be effective in imaging their structure and uncover- 
ing their true nature. We have already studied the pe- 
culiar supernova remnant G76.9+1.0 (Landecker et al. 
11391)1 and the ring G79.29+0.46 around a suspected LBV 
(Higgs et al. I1994|) . Another interesting feature is DR 16 
(Downes & Rinehart ITM)|> - G82.0+2.3 - and its sur- 
roundings. Images at 408 MHz (Paper XVIII) and 4.8 GHz 
(Wendker 1984 Paper XV) showed a cometary nebula. 
Recombination line observations at 6 cm (Piepenbrink & 
Wendker 119881 Paper XVII) gave nearly the same radial 
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velocity at three positions, suggesting that the feature is 
one physical object, with a pronounced head and an im- 
pressive flaring tail about 0.5° in length. Fig. clearly 
conveys this impression. Nevertheless, in Paper XVIII we 
labelled the feature as a ridge (CXR 10). Optically only 
the northern part of the tail is visible in Ha (filaments 
DWB 155 and DWB 157 in the catalogue of Dickel et al. 
(1969, Paper V)). Visual extinction in front of the head 
rises to more than 5 m (Dickel & Wendker ll978l Paper XI), 
totally obscuring it. 

The physical picture commonly associated with the 
term 'cometary nebula' is an object moving through its 
surroundings and producing some sort of wake. Supersonic 
motion produces a bowshock object, and is not implied in 
this case. Since sub-sonic or mildly supersonic velocities 
(less than 10 km/s) do not produce large extensions, a tail 
of length 0.5° immediately implies a very nearby object, 
or an exceptionally large one. 

In this paper we present a new VLA 20-cm image with 
resolution 15" which reveals that the apparent head of the 
comet (which we now label DR 16) has a dumbbell struc- 
ture. Furthermore, little evidence remains of a connection 
between the head and the other parts of the ridge (which 
we now refer to as the tail). We also present new observa- 
tions from the database of the Canadian Galactic Plane 
Survey (the CGPS; see Taylor et al. l2003f> and collect other 
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Fig. 2. VLA map of DR 16 at 1515 MHz. The gray scale is non-linear in order to emphasize the large-scale structure 
of the fainter emission and brightness increases from black to white. 



published observations for a comprehensive study of the 
source. In Sect. 2 to 5 we describe new observations at 
radio and infrared wavelengths and present the results. A 
distance is suggested in Sect. 6. In Sect. 7 to 9 we discuss 
the source properties and suggest a source model and a 
source history. 

2. Observations 

Inspection of Fig. ^ showed that new observations with 
resolution considerably better than 1' would be needed to 
resolve the head of the apparent comet ary nebula, with 
a substantial field of view needed to reveal the complete 



picture. In this section we will describe targeted VLA ob- 
servations and the new CGPS data products at 1420 MHz 
in continuum and 21-cm line. We also summarize results 
from other publicly available datasets. 

2,1. A VLA map at 1.5 GHz 

The observations at 1.5149 GHz were obtained with the 
NRAO Very Large Array (VLA) 1 in the C/D and D con- 
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Fig. 1. 408 MHz map of the DR 16 area extracted and 
enlarged from the map in Paper XVIII. The coordinates 
are B1950. The contours are in steps of 40 K Tb beginning 
at 140 K. 

figurations. Concurrent observations of other sources of 
interest in Cyg X have been described by Landecker et 
al. 119931 and Higgs et al. [1994 Instrumental parameters 
are summarized in Table ^ The phase centre was not 
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Table 2. Parameters of the DRAO CGPS data. 



CGPS mosaic MN2 


(release April 2002) 


Centre frequency 


1.420 GHz 


Continuum: 


Stokes I 


Synthesized beam 


71.5" x 49.1", PA -36.2° 


rms noise 


40 mK 


Line synthesized beam 


85.4" x 58.6", PA -36.0° 


rms noise 


2.5 K 


Radial velocity coverage 


-164 -» +58 km/s 


Channel spacing 


0.82446 km/s 


Velocity resolution 


1.32 km/s 



First, the head, with its size of ~5' was adequately cov- 
ered. Second, the overall structure is better depicted in 
the CGPS image discussed below. 

Fig. 121 shows a gray scale image of the whole held, 
with a grayscale chosen to emphasize the fainter large- 
scale structure. The overall impression of a cometary neb- 
ula has been lost, because the 'western ear' has become 
detached from the rest of the tail and the 'eastern ear' 
seems to be part of a long curved ridge. 

The head of the cometary nebula shows complex sub- 
structure, as seen in the detailed contour map of Fig. |3J 
Although the outer contours are nearly circular, with an 
overall size of ~5' the source has a dumbbell shape with 
some superimposed small-scale features. DR 16 is appar- 
ently placed asymmetrically on the tail, which runs north- 
south. These complex gradients make it difficult to sep- 
arate DR 16 and its components from the background 
structure. 



Table 1. Parameters of VL A observations. 



2.2. CGPS 21 cm continuum and line 



Field centre: 
a(B1950) 
<5(B1950) 
Frequency 

Date C/D configuration 
Date D configuration 
On-source time C/D 
On-source time D 
Calibrators: 

flux density 

phase 

polarization 
Synthesized beam 
Field of view (HPW) 
Synthesized map size 
Surface brightness conversion 



20 ft 30 m 45 s 
43° 37' 
1.5149 GHz 
1988 May 17 
1988 August 22 
39.7 min 
19.7 min 

3C 286; 14.39 Jy 

2050+364 

3C 138 

15" x 15" 

29.22' 

21.25' x 21.25' 

1 Jy/beam <-> 2368 K 



set on the maximum of the head but several arcmin to 
the north in order to achieve better coverage of the whole 
object. The two VLA configurations were combined and 
CLEANed using standard procedures. The map units were 
expressed in brightness temperature Tb. We used the polar 
diagram corrected map as several components appeared to 
be too extended to allow a simple post-integration factor 
to be applied. No short-spacing information was added. 



The CGPS observations include the Cyg X region. One of 
the outstanding properties of the Survey is the coverage 
of a large part of the northern Galactic hemisphere at a 
number of wavelengths with a resolution of around 1' or 
slightly worse, covering the radio and infrared regimes 2 . 
Its "backbone" is provided by the 21-cm continuum and 
line observations from the Synthesis Telescope of the 
Dominion Radio Astrophysical Observatory (DRAO). A 
description is given by Taylor (|2()()3|l . and here we present 
only those parameters relevant to the DR 16 area (in Table 
|2J). The mosaic has short-spacing information added, with 
continuum data taken from the Effelsberg survey of Reich 
et al. l|1990fl . and line data from the DRAO 26-m Telescope 
survey of Higgs & Tapping ( 2000 ). 

The CGPS continuum image of the area is presented in 
Fig. El intended to show the very complex environs of the 
alleged cometary nebula by using the scaling to empha- 
size faint emission. The area is criss-crossed by filamen- 
tary structures, which could be anything from fragments 
of shock fronts to small bubbles: none is yet identified or 
explained. The 'western ear' is now seen as a brightness 
enhancement on a filament apparently unconnected with 



2 The CGPS data are available to the astronomical commu- 
nity at |http: / /www.cadc.hia.nrc.gc.ca| 
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Fig. 3. Enlargement of the DR 16 area proper from Fig.0 The contours are in KT\>. The contours are in steps of 2 
K up to 8 K, in steps of 4 K up to 24 K, and in steps of 6 K up to 60 K. The solid, dotted or dash-dotted contours 
are assigned in order to enhance clarity. 



Table 3. Other data sets used in this paper. 



Frequency or 


Telescope 


Angular 


Reference 


wavelength 




resolution 




151 MHz 


CLFST 


102" x 70", PA 0° 


Vessey & Green [1998 


327 MHz 


WSRT 


78" x 54", PA 0° 


Rengelink et al. [TW7| 


408 MHz 


DRAO SRT 


312" x 210", PA 0° 


Paper XVIII 


4800 MHz 


MPIfR 100 m-RT 


161" x 153", PA 0° 


Paper XV 


100 \irsx 


IRAS-HIRES 


132" x 102" 


Cao et al. [TW7| 


60 /jm 


IRAS-HIRES 


102" x 60" 


Cao et al. [TWI 


25 /im 


IRAS-HIRES 


66" x 34" 


Kerton & Martin [T9U7| 


12 


IRAS-HIRES 


67" x 33" 


Kerton & Martin [15571 


8.28 


MSX 


20" x 20" 


Price 193U IRS A 
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Fig. 4. CGPS 1420 MHz grayscale image of DR 16 and surrounding region. The grayscale is non-linear, chosen to 
emphasize the large-scale structure of the fainter emission. Brightness increases from black to white. 



DR 16, and should be discussed in another context. Such 
a conclusion is not immediately possible for the 'eastern 
ear', but it seems unavoidable that we must abandon the 
picture of a cometary nebula. In this image DR 16 appears 
as a fairly compact source. Whether the 'tail' is physically 
connected to DR 16 remains an open question, to which 
we will return later. 

The CGPS 21-cm line data reveal a bewildering array 
of structure around DR 16. However, two points emerge. 
First, there are no emission structures that can be eas- 
ily related to DR 16. Second, we see that the source is 
strongly absorbed, even though its continuum brightness 



temperature is not particularly high. We will use these 
observations for a new discussion of the distance to DR 
16. 

2.3. Other observations 

Within the frame of the CGPS there are other datasets 
which yield additional information on the DR 16 area. 
We summarize these and other publicly available data in 
Table Sadly, adequate observations of the molecular gas 
are not available. The work of Cong l|1977[) and Dobashi 
et al. I)1994|l and the collection of CO survey results of 
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Leung & Thaddeus 119921 do little more than tell us that 
molecular gas exists in the general DR 16 area. Their lack 
of resolution and/or coverage prevent examination of de- 
tails. The NRAO VLA Sky Survey (NVSS) (Condon et 
al. I1998|l covers the source area. Due to its lower resolu- 
tion and insensitivity to extended structure our VLA map 
contains more fine structure. Near infrared (from 2MASS) 
and optical (POSS) images will be referenced later. 

3. Spectral energy distribution 

To establish the nature of DR 16 we began by deriving 
radio and infrared spectral indices. The difficulty was the 
complexity of the background emission. For a convenient 
and uniform approach at all frequencies we used the 'flux- 
fit' routine from the DRAO Export Software Package, 
which fits a Gaussian source component and a twisted- 
plane background simultaneously. With the present data 
this method was preferable to the alternative of defining 
a complex background and integrating the excess emis- 
sion. However convenient the technique, we do not mean 
to imply that the source consists of actual Gaussian com- 
ponents. It is clear from Figs, El that the form of the 
background and the number of Gaussians will be strongly 
resolution dependent. At a given frequency we started 
with a number of components estimated by eye, and added 
or deleted components by using the \ 2 discriminator cal- 
culated by the routine. The results of this exercise for the 
available radio data are collected in Table El The sum 

Table 4. List of sources from gaussian decomposition. 
The source notation (column 1) is frequency/component 
number. 

Source a S flux observed 



(B1950) (B1950) density size 
20 h 30 m 43° [mJy] 



151/1 


40.1 s 


29' 


33" 


1584 ± 375 


2.41' 


X 


1.36' 


151/2 


49.0 


27 


56 


1249 ± 691 


5.79' 


X 


1.20' 


327/1 


33.1 


30 


13 


759 ± 32 


2.62' 


X 


1.50' 


327/2 


33.9 


33 


01 


299 ± 31 


2.50' 


X 


1.45' 


327/3 


42.9 


31 


03 


1322 ± 41 


3.04' 


X 


2.00' 


408/1 


37.6 


30 


58 


3232 ± 327 


5.50' 


X 


4.04' 


408/2 


48.8 


33 


33 


1538 ± 468 


11.6' 


X 


3.72' 


1420/1 


31.6 


32 


16 


278 ± 18 


2.82' 


X 


1.18' 


1420/2 


32.6 


30 


19 


844 ± 19 


2.52' 


X 


1.51' 


1420/3 


42.3 


31 


09 


2105 ± 28 


3.64' 


X 


2.17' 


1515/1 


31.1 


31 


31 


78 ± 5 


1.38' 


X 


0.42' 


1515/2 


31.7 


30 


28 


1078 ± 10 


2.43' 


X 


1.28' 


1515/3 


34.5 


29 


52 


19 ± 3 


0.56' 


X 


0.31' 


1515/4 


37.3 


32 


45 


1202 ± 16 


2.92' 


X 


2.40' 


1515/5 


38.3 


30 


41 


91 ± 8 


0.67' 


X 


0.40' 


1515/6 


41.9 


31 


36 


165 ± 14 


1.34' 


X 


0.66' 


1515/7 


43.5 


30 


14 


1812 ± 17 


3.25' 


X 


2.57' 


1515/8 


44.8 


31 


02 


123 ± 13 


0.98' 


X 


0.86' 


1515/9 


47.0 


31 


32 


42 ± 5 


1.74' 


X 


0.33' 


4800/1 


38.9 


31 


03 


3805 ± 210 


4.35' 


X 


4.08' 



of the components at each frequency should be a good 



measure of the total flux density of DR 16. However, the 
spectrum constructed from these sums was rather ragged. 
We supect that the combination of angular resolution and 
frequency dependent background effects are to blame. For 
these reasons we excluded components 151/2, 408/2, and 
1515/4. The resulting total flux densities of DR 16 proper 
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Fig. 5. Radio and infrared spectral energy distribution of 
DR 16. The radio flux densities have been shifted upwards 
by a factor of 10. The broken line through the radio flux 
densities is for illustrative purposes; the high frequency 
part has a slope of —0.1. The gradient in the IR is —1.7 
and is also shown for illustrative purposes. 



are plotted in Fig. The firm conclusion from this exer- 
cise is that the radio emission of DR 16 is optically thin 
thermal emission above 1 GHz with a gradual transition 
to higher optical depth at lower frequencies. The optical 
depth of most parts exceeds unity at 151 MHz. 

The set of contour maps in the four IRAS passbands in 
Fig. reveals a general absence of structure in the source 
in the infrared. The background is relatively faint and 
smooth, and we used the 'fluxfit' routine on these data, 
just as at radio wavelengths. The resulting flux densities, 
together with one obtained in the same way from the 8.28 
/ira MSX data, are plotted in Fig. [5j 

There is obviously a non-negligible amount of cold dust 
at the location of DR 16. There is no indication that 
the IR flux density turns down at 100 /im towards the 
Planck maximum, implying a dust temperature below 25 
K, but an actual value cannot be derived; we will assign 
Idust ~ 20 K at a later point. Dust within Hn regions is 
substantially warmer (~ 40 K), and consequently the ra- 
dio and infrared parts of the spectral energy distribution 
must represent two quite different environments. 
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Fig. 6. Contour maps of the 4 infrared passbands from IRAS. Contour units in MJy/ster are: (12/im) 20, 25, 30, 
40, 50, 60, 70, 80, 90; (25//m) 30, 40, 50, 60, 80, 100, 150, 200, 250; (60^m) 200, 250, 300, 400, 500, 600, 900, 1200, 
1500; (100//m) 600, 700, 800, 1000, 1200, 1500, 1800, 2100. Dashed contours indicate a first approximation to the 
background. 



4. 21 cm continuum absorption 

The CGPS Hi data cube contains no emission features 
clearly associated with DR 16. On the other hand it is 
quite obvious that there is absorption of continuum emis- 
sion by local gas to be used for foreground and distance 
discussions. The optical depth, t, and the spin tempera- 
ture, Tg, of the absorbing neutral hydrogen column can be 
derived. For a discussion of the problems and of the equa- 
tions we refer to Wendker & Wrigge (1996, Paper XX). 

The optical depth is connected to three observed quan- 
tities via 

t{RV) = -ln(l - (T hiOB (RV) - T htOIi (RV))/T c ) 



where Tb.off, 2b j0n and T c are the expected emission line 
brightness temperature, the directly measured line tem- 
perature and the continuum brightness temperature, re- 
spectively. The variable RV is radial velocity referred to 
the local standard of rest (LSR). DR 16 is partially re- 
solved in both continuum and line images. We ignore the 
slight difference in angular resolution between the CGPS 
continuum and line mosaics (Table 

A detailed continuum map is shown in Fig. [7| where 
three on-source areas are marked. Nine-pixel averages of 
these areas are (from bottom to top): T c i = 65.0 K, 
T c2 = 59.0 K, T c3 = 56.7 K. The rms scatter of these 
averages is ~2 K, but the greatest uncertainty arises 
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DR 1 6: 1420 MHz CGPS 




82.10 82.05 82.00 81.95 

gal. longitude 

Fig. 7. Contour map of DR 16 from the 1420 MHz CGPS 
continuum data. Contours range from 14 K to 28 K in 
steps of 2 K and from 28 K to 64 K in steps of 4 K. 
The dashed squares enclose the 3x3 pixels used for 
averaging in the continuum absorption discussion (cen- 
tred at the asterisks). The outer quite regular ellipti- 
cal contours are used to define the geometric centre as 
[l,b) = (82.035°, 2.325°). 



from separation into foreground emission which is not 
absorbed and source emission which is. Other than DR 
16 itself there are at least three contributions to con- 
sider: a remnant brightness (about 5 K) from the extended 
ridge (which constituted the tail of the supposed cometary 
nebula), a diffuse Galactic component, and the cosmic 
microwave background (CMB). Although we will subse- 
quently argue that the 'tail' is not generated by the 'head', 
we nevertheless believe it is plausibly at the same distance. 
Continuum emission from the tail and from the CMB will 
therefore not confuse the absorption calculations. From 
the data of Fig. 0] we estimate that the diffuse Galactic 
component could be as large as 12 K, which could be en- 
tirely in the foreground. However, it is more reasonable to 
assume that it is spread uniformly along a 10 kpc line of 
sight through the disk. For a distance of 3 kpc to DR 16, 
the foreground contribution would be only 3.6 K. Another 
way of estimating the foreground is to associate it with the 
faint extended H n region connected with the Cyg OB2 as- 
sociation, as discussed by Huchtmeier & Wendker 1)19771 
Paper IX). Their 11-cm brightness of ^0.8 K translates 
to 3.2 K at 1420 MHz. Considering the vastly different 
assumptions, the two estimates agree quite well. We will 
subtract 3 K from the continuum brightness temperatures 
mentioned above, note that the uncertainty in this num- 
ber will subsequently dominate the absolute error budget, 
and will formally use a 4% error for T c . 



DR 16: RV = -0.64 km/s HI - CGPS 




82.15 82.10 82.05 82.00 81.95 81.90 
gal. longitude 

Fig. 8. Brightness temperature map of the three-channel 
average centred on RV = —0.64 km/s. The area is some- 
what larger than in Fig. [3 to show that the assumption 
of smooth emission on large scales is justified. A 30 K de- 
pression due to continuum absorption can be recognized 
(dashed contours). The dotted squares and asterisks have 
the same meaning as in Fig. 



The channel maps show that the scale of the H I emis- 
sion is larger than the continuum source, with emission 
features on average > 10'. A typical channel map is shown 
in Fig. IB] From inspection of this figure and Fig.0 the off- 
spectrum should be derived beyond the 100 K line contour 
and below the 30 K continuum contour. Trials with var- 
ious ring sizes showed hardly any difference between off- 
spectra obtained at different radii, and for convenience we 
chose the one closest to the source. The on-source spectra 
for the three areas and the off-source spectrum are shown 
in Fig. E| 

We now obtain optical depths and spin temperatures. 
Fig. 1101 shows 4 curves for r: the results for the individual 
areas and their mean (solid line). The error budget in the 
continuum brightness temperatures is dominated by the 
separation of foreground from background, and we esti- 
mate this to be ^4% (which may be optimistic). The un- 
certainty in on-source line temperatures comes from noise 
(receiver plus sky) and variation in emission gradients, and 
we formally set this to 5 K for a signal of ^70 K, or 7%. 
The uncertainty in the off-source spectrum is at least the 
scatter around the average, and we assign a further 7% 
to this. Conservatively adding these contributions linearly 
brings us close to 20%, which translates to an uncertainty 
in optical depth of 0.2. We therefore regard all values of 
t < 0.2 as not meaningfully measured. Looking at the 
scatter between the three individual areas we see no sig- 
nificant foreground fluctuations over DR 16 and find very 
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Fig. 9. The selected off-sepctrum is plotted as a dashed 
curve; it is the average of 72 spectra from points in a ring 
between 10 and 11 pixels away from the geometric centre. 
The three different 3x3 pixel averages (on-spectra - dash- 
dotted) are shifted downward by 20 K for clarity. Their 
average is indicated by an (unshiftcd) solid line. 




RV [km/s] 

Fig. 10. Optical depth versus the radial velocity interval 
for the local gas. The three dashed-dotted lines are for the 
individual areas and the solid line represents their mean. 

good agreement between the areas. We therefore use the 
averaged curve to derive r and T S - 

First we interpret these results in velocity space. 
Continuum absorption becomes noticeable around +19 
km/s where the emission rises to the point where T b > 30 



K. t reaches a first maximum of nearly 0.5 at +16.5 km/s, 
with Ts just below 200 K. This Hi must be the closest 
neutral layer in the Cygnus Rift, r drops slightly at +12 
km/s with an increase in T s to around 300 K, possibly 
corresponding to a neutral but warm outer layer of the 
Rift (there will be some overlap in velocity space of gas 
at distinctly different distances, and this could dilute the 
changes in r and Ts). After dropping back to 200 K, Ts 
does not vary significantly out to —5 km/s, although dif- 
ferent absorbing structures can be identified by the vari- 
ations in r. This faint anticorrelation between r and Ts 
plausibly fits into the standard picture where higher op- 
tical depth means higher volume density, which in turn 
implies lower temperature because of greater cooling ef- 
feciency. The peak in r at km/s probably marks the 
densest part of the local Rift gas. r reaches a final peak 
at —6 km/s, after which it quickly drops to unreliable lev- 
els. We intepret this as the last layer of neutral hydrogen 
in front of DR 16, which we place at —7 km/s in radial 
velocity space. 

We also computed absorption spectra for the two con- 
tinuum maxima of the tail. Absorption is detectable at 
both maxima with roughly similar optical depths, but the 
scatter rises dramatically because the continuum temper- 
ature is much lower (^35 K) and off-spectra must be de- 
rived from regions parallel to the ridge but quite far re- 
moved. Nevertheless, the results are very consistent with 
those for DR 16 proper, and from the fact that continuum 
absorption is visible out to —6 km/s we conclude that the 
tail must be at the same distance. 

The column density of neutral hydrogen can be de- 
duced by summing the relevant channels of the corrected 
emission line profile up to the continuum source. Using 
the standard relation, a channel contains A/hi = 1.823 • 
10 18 • T s {RV) ■ t(RV) ■ ARV where ARV = 0.82446 is the 
channel spacing. We obtain A/hi = 6 • 10 21 cm~ 2 . Using 
another standard relation (A/hi = 5.9 • 10 21 • £?b-v) we 
derive a total visual extinction of 3 m towards DR 16 which 
is about half the value deduced in Paper V. We suspect 
that a substantial contribution to the extinction not repre- 
sented by neutral hydrogen is connected to two molecular 
clouds which were seen in H2CO (Paper XVII) at —0.7 
and +3.4 km/s where the above relation does not hold. 
Piepenbrink & Wcndkcr (Paper XVII) estimate that the 
integrated effect of the H2CO containing clouds represents 
more than 2 m of extinction, which would close the gap. 

5. Near infrared and optical data 

5.1. MSX 

In discussing the spectral energy distribution we remarked 
that the IRAS maps do not give us much information on 
the dust distribution. The spectral index derived in that 
section implies that there is a large amount of cold dust, 
but also some warmer dust. We cannot tell whether the 
warm and cold components are mixed (usually not the 
case) or the location of the warm dust. This can be reme- 
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Fig. 11. Grayscale image of the MSX 8.28 /im map over- 
laid with radio contours from the VLA observations from 
Fig. |3| The infrared brightness increases from white to 
black. 



died to some extent by looking at the MSX data (refer- 
enced in Table |3J| . The absolute calibration of the MSX 
data is somehat uncertain, and the good fit of the MSX 
flux density with the IRAS spectrum may be a coinci- 
dence. However, the MSX image is quite useful. If we 
consider the wavelength of 8.28 /im to represent predom- 
inantly warm dust (say 40 K), then Fig. ITTI shows that, 
besides a compact concentration coinciding with the cen- 
tral radio point source, filaments of warm dust lie just 
outside the radio filaments. The picture which immedi- 
ately comes to mind is one in which the radio structures 
(ionized hydrogen) are photo-dissociation regions on the 
edges of dust clouds. We will return to this interpretation 
later. 

The MSX Archive Point-Source Catalogue contains 
nine entries over the area enclosed by the low con- 
tours of the MSX and VLA images. Seven of them co- 
incide with extended structures, and are probably not 
real point sources, although very faint point sources 
might be buried in the extended emission. The two re- 
maining entries (MSX5C G082. 0329+02. 3247 and MSX5C 
G082.0087+02.3024, later referred to as MSX8-1 and 
MSX8-2) appear to be pointlike. MSX8-2 has an optical 
counterpart in both the blue and red POSS images, and 
we regard it as a foreground star; its high brightness in 
the near-infrared photometry of the 2MASS survey argues 
for a very late type star. MSX8-1 seems to have a core- 
halo structure, and coincides with the central VLA source 
(1515/5 in Tabled}. We discuss its nature and identifica- 
tion in a later section. 



5.2. J, H, K photometry 

Two (observationally) different surveys of the DR 16 
area in the standard near-infrared bands (J, H, K) 
are available. Comeron & Torra (2001) observed at the 
Observatorio del Teide and Dutra & Bica (2001) extracted 
a dataset from the 2MASS data base. Both concluded from 
the images that a stellar over-density is present at the DR 
16 radio position and derived colour-magnitude diagrams 
(CMDs) and distances for the central parts of the pro- 
posed clusters. In view of the rather dissimilar distance 
estimates (1.8 kpc and 2.7 kpc) we obtained and examined 
both data sets. Plotting both photometry and astrometric 
positions it became apparent that the photometric values 
do agree within the mutual errors very satisfactorily. The 
main differences arc, first, a somewhat fainter catalogue 
limit for the 2MASS data, and, second, the assignment 
of a quite different area to the new cluster (in each case 
determined by eye). For the first reason we decided to use 
the 2MASS data to have an independent third look. A 
greyscale image of the cluster area is given in Fig. ^| for 
illustrative purposes. 




" RIGHT t'scEl ISIOI I !E I ?50> 

Fig. 12. K s band image of the DR 16 area from the 
2MASS data. The ellipse encloses the area which is con- 
sidered the extent of the cluster. 

We first produced a map of the stellar density distribu- 
tion in the K s band (the band least affected by extinction) 
by counting objects in 0.75 arcmin 2 circles. The central 
density enhancement was well fitted by a Gaussian, giv- 
ing a geometric centre of (B1950: 20 h 30 m 36.4 s , 43°31'01") 
and a (Gaussian) size of 1.6' X 1.1' at a position angle 
of 115°. (We do not imply that the real distribution is 
Gaussian.) The central position lies in the middle of a 
small hole, whose significance is not obvious (statistical 
fluctuation, additional patch of extinction, or division into 
two sub-structures come to mind) . The centre is about 50" 
north-west of the central VLA source (1515/5). 
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Fig. ^5] shows stellar density within elliptical rings 
around the cluster centre. The rings have a width of 0.25', 
and an axial ratio of 1.6, with the major axis at a position 
angle of 115°. The existence of a star cluster above a back- 
ground of about 8 stars per arcmin 2 is beyond doubt. The 
peak is off centre because of the hole already mentioned. 
The total extent is 3.5' along the major axis. The elliptic- 
ity may be more pronounced in the outer parts than we 
have assumed, as may be seen in Fig. El Subtracting a 
background level of 7.7 stars per arcmin 2 from the total 
of 93 stars inside this ellipse (6.01 arcmin 2 ) leaves about 
47 stars in the cluster. Only 6 stars are needed to fill the 
central depression, and the stars in the second ring appear 
to be the most heavily reddened. The central minimum is 
then probably the result of additional extinction, but the 
existence of two substructures cannot be ruled out. 



DR 16: star counts in 2MASS K-bond 




Fig. 13. Stellar density in stars per arcmin 2 plotted 
against the length of the major half-axis of elliptical rings 
with an axial ratio of 1.6:1 for the K s band. The adopted 
background level of 7.7 ±0.7 stars per arcmin 2 is the aver- 
age of all values beyond 1.75' and is indicated by a dotted 
line. The geometric centre is given in the text. 



The placement of a reddened zero-age main sequence 
on a cluster in a near infrared CMD is highly ambigious. 
The reddening vector is nearly parallel to the colour axis 
in a mK, vs. (J — K s ) CMD for the upper main sequence, 
and the distance modulus is parallel to the apparent mag- 
nitude axis. Dutra & Bica l|2001[l used a template method 
which assumes that for a small young cluster the tenth 
brightest star has an absolute magnitude close to m (spec- 
tral type about AO). The template cluster is NGC 6910 
from the optical photometry by Delgado & Alfaro ( 2000 ) . 
However, for the ninth brightest star of their member list 
we infer from the BDA database (Mermilliod 1995| a spec- 



tral type of B1V. Therefore the Dutra & Bica lj20TJTjl as- 
sumption does not seem to be well justified. 

A theoretical unreddened zero-age main sequence was 
read from the tables given by Cox J2000) for the dik vs. 
(J — K) CMD. The tables have limited coverage of spectral 
range, only from 09 to M5, but seem to be the best avail- 
able. From the parametrization of the extinction law by 
Cardelli et al. l|1989f) and values given by Mathis ( 1990J we 
calculated the reddening vector to be Ej-k = 0.168 • Ay 
where Ay is the total visual extinction to be used. The 
unreddened zero-age main sequence shifted to a distance 
modulus of 10 m is plotted in Fig.^J It will give some indi- 
cation which stars may be unreddened foreground objects 
and we see that three of the brighter stars belong to this 
group. As they are also visible on the POSS plates with 
optical colours B — V ~ 1.5 m they are indeed late type 
foreground stars. For the brighter stars there is quite a gap 
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Fig. 14. The m K vs. (J - K) CMD of all objects inside the 
elliptical area assigned to the DR 16 star cluster. The most 
leftmost line is the theoretical unreddened zero-age main 
sequence shifted to a distance modulus of 10 m , spanning 
spectral types 09V to M5V. The other lines represent the 
same sequence shifted to a distance modulus of 12. 5 m and 
successively reddened by 9 m , 12 m , 15 m and 18 m of visual 
extinction Ay. 

in colour between the unreddened sequence and the sus- 
pected cluster stars. The minimum reddening is then 1.4 m 
in J — K, corresponding to about 8 m of Ay- The scatter 
for stars with a redder colour indicates that substantial 



12 



O.P.Behre et al.: A Herbig Ae/Be star with a giant bipolar outflow in DR 16 



cluster and/or circumstellar reddening occurs. In the ab- 
sence of knowledge of the spectral type of any of the stars 
it seems impossible to 'fix' a number so that a distance can 
be inferred. The only conclusion permitted by the CMD 
is that the data are consistent with any distance between 
1 and 4 kpc. We will illustrate this with one numerical ex- 
ercise. Assume a distance modulus of 12. 5 m (3 kpc); the 
star 2MASS J2032214+434107 with m(K s ) = 9.74 m and 
(J-K s )=1.60 then has £\t-k s = 1.76 m corresponding to 
Ay — 10.5 m . Shifting it to the unreddened sequence leads 
to Mk s = — 3.9 m which in turn corresponds to a spectral 
type of 09V. Its extinction therefore occurs mostly in the 
foreground, and the additional l m of extinction could ei- 
ther be circumstellar or local to DR 16. For such a star 
the expected apparent visual magnitude is about 18. 4 m . 
But the star is visible on the red POSS plate with a mag- 
nitude around 17 m (GSC) and is just visible on the blue 
plate! And if we shift the distance modulus by l m in either 
direction we still obtain equally plausible and meaningful 
numbers. This exercise can be carried through for all stars 
above m(K s ) ~ 13. 5 m , of which the vast majority are not, 
and should not be, visible on the POSS plates, confirming 
that they are highly reddened stars, not intrinsically faint 
ones. Below this value a distinction between slightly red- 
dened foreground and highly reddened cluster stars can no 
longer be made, which explains the widely scattered band 
at the bottom of the CMD. The same arguments make it 
highly likely that at least the brighter stars in the area are 
cluster members and not field stars. Picking some sort of 
mean reddening from the CMD is not easy. 

In the end, the star cluster is only well defined by 
the projected number density distribution. Its visual fore- 
ground extinction is quite large (8 m to 9 m ) and is most 
probably produced by the Cygnus Rift at a distance of 
only 600 pc, although a distance of around 3 kpc or more 
is supported by the CMD. 

6. On the distance of DR 16 

Two different approaches are available to determine the 
distance to the DR 16 complex. On the one hand we can 
use spectral-line data for a kinematic distance. Three pub- 
lished recombination line observations are available, —7 
km/s from HUOq (Paper XVII) for DR 16 proper, -5.2 
km/s (also Paper XVII) for the tail, and —5.1 km/s from 
H87a by Lockman l)1989|l but with an angular resolution 
that enclosed DR 16 and at least part of the tail. From 
our continuum absorption study we now add —5 to —7 
km/s as a very well determined distance in radial velocity 
space. Using the rotation curve of Fich et al. (1989) this 
leads to a distance range of 2.5 to 3.5 kpc. Although kine- 
matic distances are notoriously uncertain at the Galactic 
latitudes under consideration, DR 16 seems firmly placed 
behind the Cyg OB2 association and its related features. 

On the other hand we have the two recent attempts to 
prove the existence of an open cluster connected to DR 16 
using near infrared photometry (JHK bands) and thereby 
derive a distance with the help of a colour-magnitude di- 



agram (CMD). As we showed above, the Dutra & Bica 
(2001) value of 1.8 kpc is very uncertain because the no- 
tion of a mean cluster reddening from the scatter in the 
bottom part of the CMD may be too optimistic. Comcron 
& Torra H20()l|l obtained 2.7 kpc using a different esti- 
mator (a sort of Zanstra method). However, we found 
that the star with the brightest m(K 8 ) value, 2MASS 
J2032250+434014 (see Fig. EH, was not taken into ac- 
count because it lies outside the area considered to be the 
cluster area. Adding this star into the Comeron & Torra 
( 2001 ) formalism leads to a value of 0.8 kpc, indicating the 
fragility of this method. Furthermore, the brightest cluster 
members, which contribute most of the ionizing photons, 
mostly have large circumstellar extinctions. Clearly the 
methodological problem of dust absorbing such photons 
has not been properly solved for DR 16. 

In spite of the diligent efforts in the near infrared the 
strongest constraint remains the radial velocity distance of 
—5 km/s. We maintain that this means that DR 16 is sit- 
uated behind the Cyg OB2 cluster and is physically unre- 
lated to it. We will use a distance of 3 kpc from this point, 
a choice which will appear quite plausible after we have 
discussed the nature and origin of the emission. However, 
we note that the infrared data are even better compatible 
with a larger distance. 

7. Nature of the objects 

All the data presented in preceding sections of this paper 
are compatible with the idea that DR 16 proper is an iso- 
lated cluster of very young stars. The brightest members 
can be assumed to be a group of late O to early B stars 
capable of producing enough Lyman continuum photons 
to ionize the observed Hn. However, the large number 
of stars with substantial circumstellar extinction creates a 
difficulty. Why are no ultracompact H n regions (UCHIIR) 
detected? If any exist they must be so compact that they 
can be observed only at higher frequencies. This would 
then rule out any substantial leakage of Lyman continuum 
photons. Furthermore, UCHIIR should be seen as point 
sources at infrared wavelengths (e.g. Wood & Churchwell 
1989), which is not the case. The dust distribution seen 
in the infrared seems to be rather diffuse, except for the 
warm dust, which may be concentrated in filaments on the 
outskirts. 

The ease with which the VLA map can be decomposed 
into Gaussian components could indicate that there are 
multiple ionization centres. Indeed, each small component 
needs only a middle B type star and the large diffuse emis- 
sion a single 09V one. However, we would then expect to 
see these stars more or less centred on the small radio 
components. The fact that only one such identification 
appears from a careful overlay, namely VLA 1515/5 with 
2MASS J2032213+434056 (which we refer to together as 
#5) argues against this geometry. We suggest that the 
radio structures seen are density fluctuations, with the 
exception of #5. 
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A first estimate of the number of Lyman continuum 
photons required for the ionization of the Hn region DR 
16 can be made using 

N Lc = 4.761 • 10 48 • T -°- 45 ■ zA 1 ■ S(v) ■ D 2 s" 1 . 
This relation (Mezger I1972fl is valid only for an optical 
depth less than unity, which appears to be the case. Using 
an electron temperature of 10 K, an average frequency 
of 1.45 GHz, a total observed flux density of 3.3 Jy, and a 
distance of 3 kpc we obtain ATl c = 2.3 • 10 48 s _1 . Since we 
will subsequently suggest that the three Gaussian compo- 
nents, #5 (the central VLA source), 1515/3, and 1515/6, 
which define the central dumbell appearance of DR 16, are 
a Herbig Ae/Be (HAeBe) star with a giant bipolar out- 
flow, this value represents an upper limit. We now inspect 
the brightest cluster stars to see whether they can deliver 
so many Lyman continuum photons. The two stars bright- 
est in apparent m K magnitude (2MASS J2032250+434014 
and 2MASS J2032126+434103) are also heavily reddened 
by circumstellar material (Fig. I14J) . This additional ex- 
tinction amounts to 8.5 m and 13 m , respectively. Although 
their position in the CMD makes them late O-type stars, 
they will not contribute to ionization of the diffuse H II, 
since we consider that they are still within their parent 
molecular cocoons. The absence of UCHIIR around them 
supports this suggestion - the UCHIIR are still too small 
and our frequency is too low. The third brightest clus- 
ter member is 2MASS J2032214+434107, already used 
as an example in Sect. 5 and shown to be most likely a 
B0.5V star. Its Lyman continuum photon production is 
8 • 10 47 s _1 (Vacca et al. 1199 6). This star is marginally ad- 
equate, and, if supplemented by the next ranked B stars 
(type B2), would account for the ionization of the diffuse 
Hn. The star lies near the centre of the object, making it 
plausible that DR 16 is an ionization bounded H n region. 
Remarkably, increasing the distance modulus to 13 m (4 
kpc) inverts problem, providing too many photons (which 
of course could be absorbed by dust). 

The two next brightest stars (2MASS 
J2032214+434107 and 2MASS J2032197+434129) 
should have 4 m and 6 m of circumstellar (or local) extinc- 
tion according to the CMD and are consequently invisible 
or barely visible on the POSS. 

The seventh brightest (2MASS J2032213+434056) star 
deserves special attention. It is the only one which co- 
incides with one of the radio components, the small 
Gaussian component 1515/5 fitted to the central point 
source in the VLA map. It is also the only one to co- 
incide with a near-infrared source (see Figs. 03 and 1111 
#5). Working through the same exercise as above, (mK = 
11. 7 m , A w = 10.4 m -> M K = -1.9 m -> type B2V 
My = —2.4) we calculate an expected apparent red mag- 
nitude of 17. 8 m . This appears to be compatible with the 
GSC red brightness of 16. 7 m , given the roughness of the 
calculations. We conclude that #5 has not yet reached 
the main sequence and thus qualifies as a HAeBe star. In 
addition we interpret the coincidence with a small-scale 
radio source and the dumbbell shaped structure in the ra- 
dio map as an indication of a central source plus jet. The 



radio lobes are then remnants of a giant bipolar outflow. 
Lower mass stars in a less photoionizing surrounding often 
produce Herbig-Haro objects (HH) in such circumstances. 
The projected outflow length on either side is about 1' 
or not quite 1 pc. In this interpretation we assume that 
the diffuse Hn and the dumbbell are projected onto each 
other without physical interaction at the present time, but 
both are features of the star formation event of DR 16. 

It is one thing to make this qualitative interpretation, 
but to make it more quantitative is quite another. In the 
following discussion we take properties and relations for 
HAeBe stars from the review of Waters & Waelkens (Q.998 ) 
and for HH objects from that of Reipurth & Bally l|2001|l . 

The requirements for classification as a HAeBe star are 
only partially met. Although it seems certain that this is 
a B star, we lack a suitable spectrum to test the emis- 
sion line requirement; we need a high-resolution spectrum 
of a fairly faint star. On the other hand, the infrared ex- 
cess is definitely there. The existence of a point source 
in the MSX 8.3 /Ltm band and its visibility at other MSX 
wavelengths (Fig. 111(1 is a strong indicator. Both the ra- 
dio source and its near-infrared counterpart are definitely 
elongated in the direction of the alleged jet structure but 
hardly perpendicular to it. This is best measured in the 
radio domain by the sizes of the fitted Gaussian: the ma- 
jor axis is 40" and the minor axis less than 20". However, 
some sort of core-halo geometry is equally compatible with 
the observations. We emphasize that the data do not have 
enough angular resolution to separate emission from the 
(extended) atmosphere of the star, from the remaining na- 
tal disk, or from the polar jet. Similarly, we cannot really 
say whether a jet connects the central source and the bipo- 
lar plumes. Bipolar outflows are often detected from their 
CO emission, and we have already lamented the lack of 
adequate observations of this kind. We suggest that the 
radiation from the emergent star has mostly ionized the 
outflow plumes. The warm dust, visible at 8.3 /jm, is po- 
sitioned at the outside of the Hn as seen from the star 
(vividly portrayed in Fig. Illfl and we suggest that this 
large surface is the working interface between the outflow 
and the remaining cold interstellar medium from the par- 
ent molecular cloud. 

8. Remarks on the 'tail' 

H I absorption of the continuum emission from the tail is 
detected albeit very weak out to the same radial veloc- 
ity as for DR 16 (see Sect. 4). Recombination lines were 
found at similar velocities (Sect. 6). Landecker 1)1984[) de- 
tected a H166a line with a 0.6° beam indicating that 
widespread H n with slightly negative radial valocities 
must be present. Cong H1977[) found a molecular cloud 
extending 44' by 20' with a radial velocity of about —2.5 
km/s. All these data provide circumstantial evidence that 
the tail should lie at about the same distance as DR 16, 
but none imply any physical interaction. 

In the original low-resolution image (Fig.Q the source 
flares out considerably. It is apparent from Fig. 0] that 
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this impression is created by long filaments with varying 
brightness which do not cross DR 16 implying that they 
are not part of it. The western arc is narrow with a radius 
of at least 0.5° ; its outline is lost in the general confusion 
in the area. A working hypothesis is that a shock of distant 
origin ran into the denser medium of the DR 16 complex, 
producing the filamentary structure. For the tail itself, we 
suggest that it is part of a ring-like feature with a radius 
of about 6' (Figs. Eland 01, well defined by emission from 
half of its perimeter, and by an apparent central hole with 
a linear extend of about 10 pc. 

More striking than the ring itself is the point source at 
a (1950) = 20 h 31 m 40.82 s , <5i 950 = 43°43'34.7", very close 
to its centre. However, the point source seems to be non- 
thermal, with flux densities of 117 ± 1 mJy at 1515 MHz, 
113 ± 6 mJy at 1420 MHz, 249 ± 15 mJy at 327 MHz 
and 879 ± 202 mJy at 151 MHz, indicating a spectral 
index of around —0.8. This would be compatible with an 
extragalactic source such as an AGN or a quasar. 

While examining the optical appearance of the area 
we noticed that the three stars brightest in mx are all 
about 1' south of the centre of the ring (at B1950 positions 
of 20 h 31 m 40.85 s /43 o 42'22.0", 20 h 31 m 45.31 s /43°42'11.7", 
and 20 h 31 m 39.35 s /43°41'43.1" ). Their colours indicate 
heavy reddening. From the magnitudes given in the 
USNO-B and 2MASS catalogues they could be 05 to Bl 
stars at 3 kpc distance with a reddening of > 6 m , just 
the amount of foreground reddening for the DR 16 clus- 
ter. The ring may have been created by the action of the 
stellar winds of these stars on their surroundings, but we 
have not proved any connection. 

9. Summary and discussion 

We began this study with the idea that we were seeing an 
object which seemed to be a very large cometary nebula. 
In summary the higher resolution data have revealed that 
we were looking at a complex of different objects. The 
main radio continuum source, now identified as DR 16 
proper, is the common Hn region of an open star cluster. 
The cluster is small, with a total mass of a few hundred 
solar masses, and a substantial fraction of the low-mass 
stars have yet to reach the main sequence. On the other 
hand, the fact that at least some of the O stars have shed 
their natal envelopes and are ionizing a common Hil en- 
velope points to an age of some hundred thousand years. 
However, DR 16 is not an isolated group. We have found 
several signs of interaction with surrounding material. One 
of the neighbouring features, in our original picture the 
'tail' of the cometary nebula, lies at roughly the same dis- 
tance and is probably interacting with the same molecular 
cloud as DR 16. 

We used various methods to estimate the distance to 
the complex and derived a best value of 3 kpc. The uncer- 
tainty on this number is probably 0.5 kpc, with the near 
distance constrained by the properties of the stellar clus- 
ter and the far one by the radial velocity of the H n region. 



This distance is large enough that physical interaction be- 
tween DR 16 and Cyg OB2 seems quite implausible. 

Despite having many new observations, the lack of ad- 
equate CO data impedes interpretation. Cong l|1977|l es- 
timated that there is about 85 000 Mq of molecular ma- 
terial in the area (scaled by us for distance). What is re- 
quired is a fully sampled CO map with a resolution of 1' 
or better. We predict that such a map would show much 
small-scale structure. 

The internal structure of the DR 16 H II region is rather 
smooth, with only minor density fluctuations after the 
components 1515/5, 1515/3 and 1515/6 are subtracted. 
These three components together define the structure we 
have called the dumbbell. We have shown that enough ion- 
izing photons are available from stars on or near the main 
sequence to produce a classical Hn region. The lack of 
UCHIIR around those stars with substantial circumstel- 
lar material could be an age effect. This idea is supported 
by the recent considerations of Keto H2002|) who has ar- 
gued that in the beginning the ionization front will stall 
at a very short distance from the star for a rather long 
time. 

This leaves the dumbbell as the most interesting in- 
dividual object. Furthermore, it has an early type star 
within its central radio source. That the dumbbell is a 
single physical object is not necessarily obvious. The best 
evidence is perhaps in Fig. II II where the two outer compo- 
nents appear to be confined by the only two well-defined 
substructures in the dust distribution. Furthermore, the 
central radio source is definitely elongated in the direc- 
tion of the lobes. Together these clues suggest a bipolar 
outflow (keep in mind that the torus required for the flow 
is normal to the dumbbell axis). 

10. Conclusions 

Part of the following conclusions are based on two review 
papers, one on UCHIIR by Churchwell (2002), and one 
on Herbig Ae/Be stars by Waters & Waelkens (1998). 
Massive stars produce UCHIIR which are dispersed after 
the stars reach the main sequence. The oldest, most mas- 
sive stars in the DR 16 cluster have nearly all reached this 
state and can now provide the ionization of the dispersed 
diffuse gas. The younger ones are far enough from the 
main sequence that their small UCHIIR are not observ- 
able. As the lines corresponding to stellar evolution and to 
decreasing circumstellar extinction are nearly parallel in 
the CMD we have available (see Fig. 114ft - we cannot pro- 
vide any more details of this picture. The lower mass limit 
for young stars to which this picture applies is usually put 
around 10 Mq (or spectral type B0. 5). Less massive pre- 
main-sequence stars are usually called HAeBe stars and 
develop more slowly towards the main sequence, but can 
shine through their natal cloud much sooner since their ac- 
cretion disks are less massive. Probably as a consequence 
of this, outflows are less vigorous and hardly observed at 
stages near the main sequence. Images of the thermal ra- 
dio emission portray the average electron density along the 
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line of sight, and this can smooth out small-scale struc- 
tures. Remains of old UCHIIR and of bipolar outflows, 
which are probably ubiquitous in such a region, are thus 
not visible in our radio maps. 

Our suggestion that the central star of the source 
1515/5 is a HAeBe star and that the Hn structures 1515/3 
and 1515/6 are the working surfaces of the lobes of a bipo- 
lar outflow impinging on neutral and/or molecular mate- 
rial implies the following specific evolutionary status. The 
star is close enough to the main sequence that typical 
HAeBe properties are weak (this statement can be tested 
by spectroscopy) . Its mass is large enough that it also has 
properties usually ascribed to the massive star picture i.e. 
an earlier strong outflow which has now ceased. On the 
other hand, the outflow is sufficiently recent that the den- 
sity fluctuations it produced are still visible, although its 
ionization may by now be maintained by the diffuse radi- 
ation rather than by the originating star. This transition 
status makes the star quite interesting. 

On the grander scale, DR 16 can now be added to the 
list of sources in Cyg X which comprise an Hn region 
and a small stellar cluster and lie at about a distance of 3 
kpc. This is beyond the zone of influence of the Cyg OB2 
association, often regarded as all dominating. The Cyg 
OB2 association was suggested to be a young globular 
cluster rather than a classical open cluster (Knodlseder 
2000). Distances around 3 kpc have been suggested for 
DR 21 and DR 23 (Papers XVII and XI, respectively). 
The projected distance of DR 16 from these is only 100 pc 
or so, which may indicate that this part of the local spiral 
arm contains an extended complex of small star forming 
regions. DR 16 itself is quite bright, about one third of 
the Orion nebula. Compared to DR 15, itself a bright and 
conspicuous source in Cyg X, it is about 15 times brighter. 
On the other hand, it has only about l/7th of the radio 
power of DR 21, usually regarded as a prominent cluster 
of just forming stars. 
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